close all; clear all; path(pathdef); clc;

% add subfunctions
addpath(genpath('subfunctions/'))
addpath(genpath('COMSOL_simulation/'))
addpath(genpath('data_suppinfo/'))

% color library
color_list{1} = [67 84 147]/255;        % dark blue
color_list{2} = [110 153 201]/255;      % light blue
color_list{3} = [196 40 27]/255;        % dark red
color_list{4} = [218 134 121]/255;      % light red
color_list{5} = [40 127 70]/255;        % dark green
color_list{6} = [161 196 139]/255;      % light green
color_list{7} = [218 165 32]/255;       % gold
color_list{8} = [238 221 130]/255;      % light gold
color_list{9} = [128 0 128]/255;        % purple
color_list{10} = [204 153 255]/255;     % light purple
color_list{11} = [161 165 162]/255;     % grey

%% Load data
load('glass.mat');
load('embryo.mat');

bPlotOption.MarkerSize  = 4;
bPlotOption.CapSize     = 0;
bPlotOption.LineWidth   = 1;
bPlotOption.FontSize    = 12;
bPlotOption.MarkerList = {'o', 's', 'd', '>', 'x'};

fh = figure(1); cla()
set(fh, 'position', [100 100 [180 230]])
cla()

% simulation
plot(glass.x_sim,glass.T_sim, 'color', color_list{5}, 'linewidth', 1)
hold on
plotcustom(glass.x_exp, glass.y_exp, glass.yerr_exp, bPlotOption, 1, 5)

plot(embryo.x_sim,embryo.T_sim, 'color', color_list{9}, 'linewidth', 1)
plotcustom(embryo.x_exp, embryo.y_exp, embryo.yerr_exp, bPlotOption, 2, 9)

axisinfo('Distance (\mum)', '\DeltaT (K)', [0], [0], bPlotOption.FontSize)
set(gca, 'YTick', [0 10 20])